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摘要 :【 目的 】 本 研究 旨 在 通过 趋势 分 析 对 胁迫 意大利 蜜蜂 Apis mellifera ligustica ( 简称 意 蜂 ) 4i 3 
Vi xd Ascosphaera apis 的 差异 表达 基因 (DEGs) 进 行 转 录 组 分 析 。【 方 法 】 将 纯化 的 球赛 菌 
孢子 配制 为 1 x10" 孢子 [mL 的 饲料 饲 喂 意 蜂 3 日 龄 幼虫 ,利用 Jumina HiSeq 2500 P&R d 
胁迫 的 意 蜂 幼虫 肠 道 cDNA. 进行 测序 ,将 过 滤 得 到 的 有 效 读 段 (clean reads) 映射 (mapping) 到 核糖 
体 数据 库 及 意 蜂 参考 基因 组 ,最 后 将 未 映射 上 的 reads 映射 到 本 课题 组 组 装 并 注释 的 球 吉 菌 参考 
转录 组 。 利 用 STEM 软件 对 DEGs 进行 趋势 分 析 。 利 用 WEGO 软件 对 显著 表达 模式 中 的 DEGs 3t 
ff CO 富 集 分 析 。 利 用 Blastall 对 显著 表达 模式 中 的 DEGs 进行 KEGG 代谢 通路 富 集 分 析 。 最 后 ， 
通过 对 随机 选取 的 6 个 DEGs 进行 RT-qPCR 分 析 , 对 RNA-seq 数据 进行 验证 。【 结果 】 球 澡 菌 胁迫 
意 蜂 幼虫 肠 道 样品 的 Jlumina W) EP 2-43: 21] 3 d A 83 25 454 076 条 原始 读 段 (raw reads) ,经 过 滤 得 
到 24 909 820 条 clean reads,Q30 34 Æ 93. 4696 以 上 。 趋 势 分 析 结 果 显 示 ,19 893 个 DEGs 聚 类 为 8 
个 表达 模式 ,其 中 有 12 151 个 DEGs 聚 类 为 3 个 表现 为 显著 上 调 趋势 的 表达 模式 。G0O 富 集 分 析 结 
果 显 示 , 表 现 上 调 趋势 的 DEGs 富 集 在 40 个 GO term, 富 集 基 因数 最 多 的 为 细胞 进程 (cellular 
process) (2 601 unigenes) ,其 次 为 代谢 进程 (metabolic process) (2 553 unigenes) 和 细胞 (cell) (2 522 
unigenes) o KEGG 代谢 通路 富 集 分 析 结 果 显 示 , 上 调 趋势 中 的 DEGs 富 集 在 119 个 代谢 通路 上 ,其 
中 富 集 基因 数 最 多 的 是 核糖 体 (ribosome ) (213 unigenes) ,其 次 为 氨基 酸 生 物 合成 (biosynthesis of 
amino acids) (154 unigenes) 和 内 质 网 蛋白 加 工 (protein processing in endoplasmic reticulum ) ( 130 
unigenes)。 共 有 48 个 DEGs 富 集 在 MAPK 信号 通路 上 , 聚 类 分 析 结 果 显 示 , 这 些 DEGs 随 着 胁迫 
时 间 的 延长 表达 水 平 逐渐 增强 。RT-qPCR 结果 显示 ,6 个 DEGs 的 表达 水 平 变化 趋势 与 RNA-seq 
数据 一 致 ,证 明了 本 研究 中 的 转录 组 数据 真实 可 靠 。【 结论 】 本 研究 为 在 分 子 水 平 揭示 球 才 菌 的 致 
病 机 理 提 供 了 重要 信息 ,也 为 阐释 球 吉 菌 胁迫 意 蜂 幼虫 过 程 中 的 病原 -和 宿主 互 作 机 制 芮 定 了 基础 。 
关键 词 : 意大利 蜜蜂 ; REA; 幼虫 肠 道 ; RNA-seq; 转录 组 ; 差异 表达 基因 

中 图 分 类 号 : Q966 ”文献 标识 码 : A ”文章 编号 : 0454-6296(2017)04-0401-11 
Transcriptomic analysis of Ascosphaera apis stressing larval gut of Apis 
mellifera ligustica ( Hyemenoptera: Apidae) 

CHEN Da-Fu'*, GUO Rui ** , XIONG Cui-Ling"*, LIANG Qin', ZHENG Yan-Zhen', XU Xi-Jian!, 
HUANG Zhi-Jian', ZHANG Zhao-Nan', ZHANG Lu', LI Wen-Dong' , TONG Xin-Yu', XI Wei-Jun? 
(1. College of Bee Science, Fujian Agriculture and Forestry University, Fuzhou 350002, China; 
2. Jinhua Academy of Agricultural Sciences, Jinhua, Zhejiang 321000, China) 















































& 


基金 项 目 : 国家 自然 科学 基金 项 目 (30800806 ) ; 现代 农业 产业 技术 体系 建设 专项 资金 (CARS-45-KXJ7 ) ; 福建 农林 大 学 科技 发 展 资 4 
(KF2015123) 

作者 简介 : 陈 大 福 , 98, 1973 年 3 月 生 , 河南 信阳 人 , 副教授 , 研究 方向 为 蜜蜂 保护 学 , E-mail; dfehen826 163. com; E, 9, 1987 年 5 月 
Ak, 安徽 六 安 人 , 讲师 , 研究 方向 为 蜜蜂 保护 学 , E-mail: fafu_ruiguo@ 126. com; ERI, 女 , 1977 年 1 HÆ, 福建 龙岩 人 , 实验 
Jf, 研究 方向 为 蜜蜂 保护 学 ,E-mail: 992355037@ qq. com 

# 共 同 第 一 作者 Authors with equal contribution 

“通讯 作者 Corresponding author, E-mail; fafu_ruiguo@ 126. com 

收 稿 日 期 Received; 2016-12-26 ; 接受 日 期 Accepted; 2017-02-23 












































402 


昆虫 学 报 Acta. Entomologica Sinica 


60 2$ 





Abstract: [ ^im] This research is designed to conduct transcriptomic analysis of the differentially 
expressed genes ( DEGs) of Ascosphaera apis stressing larval guts of Apis mellifera ligustica via trend 
analysis. [Methods] The purified A. apis spores at a concentration of 1 x 10” spores/mL was used to feed 
3-day-old larvae of A. m. ligustica, and then the cDNA of stressed larval guts was sequenced at Illumina 
HiSeq 2500 platform. After filtration, the clean reads were used to mapping the ribosome database and 
the reference genome of A. m. ligustica, and the unmapped reads were used to mapping the reference 
transcriptome of A. apis assembled previously. The STEM software was used to analyze the gene 
expression patterns. Gene ontology ( GO) enrichment analysis for DEGs involved in significant expression 
patterns was performed using WEGO software. KEGG enrichment analysis for DEGs associated with 
significant expression patterns was carried out by using Blastall. Finally, RT-qPCR analysis of six 
randomly selected DEGs was performed to validate the RNA-seq data. [Results] The RNA-seq of A. apis 
produced 25 454 076 raw reads, and after filtration, 24 909 820 clean reads with Q30 above 93. 4696 
were obtained. Trend analysis results showed that 19 893 DEGs were grouped into eight gene expression 
patterns, among which 12 151 DEGs were assigned to three expression patterns with significantly up- 
regulated expression trend. GO enrichment analysis results indicated that all DEGs within significantly 
up-regulated expression patterns were enriched in 40 GO terms, and the mostly enriched one was cellular 
process (2 601 unigenes) , followed by metabolic process (2 553 unigenes) and cell (2 522 unigenes). 
KEGG enrichment analysis result displayed that all DEGs within significantly up-regulated expression 
patterns were enriched in 119 metabolism pathways, and the mostly enriched one was ribosome (213 
unigenes) , followed by biosynthesis of amino acids (154 unigenes) and protein processing in 
endoplasmic reticulum (130 unigenes). Furthermore, it was found that 48 DEGs were enriched in MAPK 
signaling pathway, and the cluster result suggested that the expression levels of these DEGs increased as 
the stressing time prolonged. RT-qPCR results demonstrated that the expression patterns of the six DEGs 
were consistent with those of RNA-seq data, confirming that our transcriptome data are credible. 
[Conclusion] The findings in this study not only provide the key information for uncovering the 
pathogenesis of A. apis at the molecular level, but also lay a foundation for clarifying the pathogen-host 
interaction in A. m. ligustica under the stress of A. apis. 

Key words: Apis mellifera ligustica; Ascosphaera apis; larval gut, RNA-seq; transcriptome; differentially 
expressed genes ( DEGs) 


据 报 道 ,蜜蜂 为 全 球 7096 的 作物 和 野生 植物 授 
粉 (Klein et al., 2007; Elke, 2010)。 蜜 蜂 易 遭受 多 
种 病原 危害 而 雏 患 细菌 病 、 真 菌 病 、 病 毒 病 等 疾病 
2010), 其 中 ,白垩 病 
(chalkbrood) 是 一 种 最 具 代 表 性 的 蜜蜂 真菌 病 ,1913 
年 Massen 在 德国 首次 报道 了 对 白垩 病 的 观察 
( Aronstein and Murray, 2010) ,中 国 大 陆 1990 年 发 
生 白 亚 病 ( 梁 勤 和 陈 大 福 ,2008 ) 。 近 年 来 , 随 着 蜂 
产品 贸易 的 全 球 化 ,白垩 病 的 发 病 率 逐年 上 升 
(Aizen et al., 2009), H 3E jg Ji p Z E EK 3E TDI 
Ascosphaera apis 特异 性 侵 染 蜜蜂 幼虫 而 导致 的 疾 
病 ,发 病 时 期 一 般 为 春季 和 初夏 , 虽 不 对 蜂 群 造成 毁 
灭 性 打击 ,但 却 能 造成 成 年 蜜蜂 数量 的 大 幅 下 降 , 从 
而 严重 影响 蜂蜜 等 产品 的 产量 ( Bailey, 1963 ) , 据 报 
道 ,该 病 可 造成 蜂蜜 产量 下 降 5% ~37% ( Zaghloul 
et al., 2005) 。 国 内 养 蜂 生产 中 使 用 的 主要 蜂 种 是 























( Evans and Spivak, 



























































意大利 蜜蜂 Apis mellifera ligustica ( 属于 西方 蜜蜂 A. 
mellifera , 简称 意 蜂 ) 和 中 华 蜜蜂 Apis cerana cerana 
(属于 东方 蜜蜂 4，cerana, 简称 中 蜂 ) ,前 者 极 易 被 
PRE VR [3 TUE TEE A. EA ES , 而 后 者 具有 和 较 强 的 球 嘻 
菌 抗 性 。 球 圳 菌 对 蜜蜂 幼虫 的 侵 染 主要 是 通过 蜜蜂 
的 取 食 活动 , 含 球 者 菌 孢 子 的 食物 被 蜜蜂 摄 和 人体 内， 
孢子 在 蜜蜂 消化 道 萌 发 生长 ,继而 突破 消化 道 围 食 
膜 和 中 肠 屏 障 ,在 体腔 大 量 生长 繁殖 ,导致 蜜蜂 幼虫 
死亡 ,最 后 突破 体 表 , 异 宗 菌 丝 接触 形成 大 量子 代 了 
子 (李江 红 等 , 2012)。 

yr 20 年 来 ,国内 外 学 者 在 培养 方法 (Anderson 
and Gibson, 1998; Laf，2015 ) 形态 学 (Spiltoir， 
1955; Skou, 1988; Anderson and Gibson, 1998) . 病 
理学 ( Bailey and Ball, 1991; 
Chantawannakul, 2008), 、 流行 病 学 (Flores et al., 
2005a，2005b )、 增 殖 方 式 ( Póggeler, 2001; 
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Chorbifiski, 2003) 、 蜜 蜂 防 御 (Stanley et al., 2009; 
Evans and Spivak, 2010) 以 及 疾病 防治 (Hornitzky , 
2001; Tarpy, 2003; Mourad et al., 2005 ) 等 方面 对 
HERFST RIAR AAH EERE A Y 
病理 和 检测 方面 开展 了 较为 系统 的 研究 , 如 梁 勤 等 
(2001) MERKIR AI EER .矿质 元 素 等 方面 研究 
了 营养 生态 条 件 对 球面 菌 生长 及 产 孢 的 影响 ,结果 
表明 营养 生态 条 件 的 变化 对 球 早 菌 的 影响 极 大 ; 郑 
志 阳 等 (2011 ) 对 健康 和 患 病 蜜蜂 幼虫 血 淋 巴 进行 
SDS-PAGE 电泳 和 蛋白酶 、 酯 酶 的 活性 染色 ,发 现 健 
康 蜜蜂 幼虫 血 淋 巴 中 的 蛋白 含量 丰富 , 主要 由 4 种 
高 分 子 质 量 的 蛋白 组 成 ,而 患 病 幼虫 血 淋 巴 中 的 蛋 
白 含量 很 少 , 主要 蛋白 组 分 被 降解 ,多 种 蛋白 酶 和 酯 
酶 的 活性 在 患 病 幼 虫 血 淋巴 中 可 以 检测 到 , 而 在 健 
康 幼虫 中 检测 不 到 。 然 而 ,由 于 基因 组 信息 的 缺失 ， 
球 喜 菌 在 分 子 水 平 的 研究 进展 缓慢 。2006 年 Qin 
TAM T EREA ARSEF7406 菌株 的 基因 组 序列 
(Qin et al., 2006) , JJ ERSETRLRS 2) T AARE T 3E 
础 ,但 作者 当时 并 未 公布 基因 的 位 置 和 功能 注释 信 
息 。Cornman 等 (2012 ) 利 用 Roche 454 焦 磷 酸 测序 
技术 对 来 自 培养 基 的 球面 菌 菌 丝 和 来 自 僻 蜂 幼 虫 感 
染 组 织 的 球 时 菌 菌 丝 进行 了 转录 组 测序 ,功能 分 析 
表明 球 圳 菌 的 差异 表达 基因 (differentially expressed 
genes, DEGs) 参与 了 交配 类 型 .细胞 内 信和 号 转 导 和 
应 激 反应 。 

目前 ,有 关 球 喜 菌 的 组 学 研究 极 少 ,白垩 病 过 程 
中 的 病原 -宿主 互 作 机 制 尚 不 明确 。 前 期 研究 中 ,我 
们 已 成 功 组 装 球 喜 菌 的 参考 转录 组 并 对 其 进行 功能 
及 代谢 通路 注释 ( 张 盟 楠 等 , 2017 ) ,可 为 利用 胁迫 
意 蜂 幼虫 肠 道 的 球 宫 菌 的 转录 组 学 研究 提供 重要 的 
参考 信息 。 本 研究 在 前 期 研究 基础 上 对 胁迫 意 蜂 幼 
虫 肠 道 的 球 早 菌 的 差异 表达 基因 (DEGs) 进行 转录 
组 分 析 ,研究 结果 可 为 阐释 球 宫 菌 的 致 病 机 理 及 球 

























































































德国 Qiagen 公司 , 高 碘 酸 钠 购 自 美 国 Sigma 公司 ， 
DNA Ligase 购 自 美国 Thermo 公司 ,RNAiso Reagent 
jid AG] x. Ex Taq Polymerase 及 Superscript M 
Reverse Transcriptase 均 购 自 日 本 TaKaRa 公司 。 
cDNA 文库 构建 试剂 盒 TruSeq™ DNA Sample Prep 
Kit -Set A 为 美国 Jumina 公司 产品 。 其 他 试剂 均 
为 国产 分 析 纯 。 

恒温 恒 湿 气候 箱 购 自 中 国 宁 波江 南 仪器 厂 , 超 
净 工 作 台 购 自 中 国 苏州 安泰 空气 技术 有 限 公 司 , 超 
低温 冰箱 购 自 合肥 中 科 美 将 低温 科技 股份 有 限 公 
F] ,实时 效 光 定量 PCR 仪 (ABI 7500) 为 美国 Life 公 























司 产品 。 
1.3 意 蜂 幼 虫 饲养 与 球 圳 菌 孢子 纯化 





意 蜂 幼虫 的 饲料 配方 参照 王 倩 等 (2009 ) 的 方 
法 配制 。 预 实验 结果 显示 7 日 龄 意 蜂 幼虫 的 成 活 
率 可 达 70% 以 上 。 将 2 日 龄 幼虫 移 至 24 孔 板 培 
养 板 ,35Y ,相对 湿度 70% 条 件 下 饲养 ,每 24 h 更 
换 饲 料 。 

将 4 保存 的 球 早 菌 培养 迄 置 于 已 紫外 灭 菌 的 
超 净 台 , 用 酒精 灯 上 为 烧 片 刻 的 锯 子 夹 取 少量 菌 丝 
在 平板 中 央 2 cm 左右 圆心 区 域内 划 线 操作 , 盖 上 培 
养 耻 , 置 于 37% 生化 箱 恒 温 培养 ,3 d 后 观察 球 圳 菌 
生长 情况 。 接 种 8 ~ 10 d 以 后 , 待 培养 由 上 黑色 和 孢 
子 较 多 时 ,将 孢子 刮 至 干净 的 EP 管 中 , 充 分 研磨 后 
按照 Jensen 等 (2013 ) 离心 纯化 ,梯度 稀释 后 用 血球 
计数 板 对 孢子 进行 计数 。 
1.4 深度 测序 样品 准备 

将 高 浓度 的 球 洁 菌 纯化 孢子 混入 人 工 饲 料 , 稀 
释 至 1 x 10" 孢子 /mL, 饲 喂 3 日 龄 幼虫 ,4 日 龄 开始 
饲 喂 不 含 球 洁 菌 孢 子 的 正常 饲料 。 分 别 剖 取 球 时 菌 
胁迫 的 4, 5 和 6 日 龄 意 蜂 幼虫 肠 道 样品 (分 别 记 为 
AamT1 AamT2 和 AamT3 ) , 液 氮 速冻 后 迅速 转移 至 
超低温 冰箱 保存 备用 。 上 述 各 样品 均 设 3 个 生物 学 


















































赛 菌 - 意 蜂 幼虫 的 互 作 机 制 提供 重要 信息 ,也 为 白垩 
病 的 防治 莫 定 基础 。 


1 材料 与 方法 


1.1 供 试 蜜蜂 幼虫 与 球 圳 菌 

本 研究 使 用 的 意 蜂 幼虫 取 自 福建 农林 大 学 蜂 
学 院 教学 蜂 场 , 球 者 菌 菌 株 由 福建 农林 大 学 蜂 学 
院 蜜蜂 保护 实验 室 保 存 并 活化 。 
1.2 主要 实验 试剂 及 仪器 

DNasel 和 Oligotex mRNA Kits Midi 试剂 盒 购 自 
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重复 。 
1.5 cDNA 文库 构建 及 了 lumina 测序 

利用 RNAiso Reagent 试剂 盒 抽 提 1.4 节 各 样品 
的 总 RNA ,然后 用 RNase-free DNase 工 去 除 基 因 组 
DNA 残留 。RNA 的 质量 通过 琉 脂 糖 凝 胶 电 泳 和 
NanoDrop ND-2000 ( NanoDrop, Wilmington, DE, 
USA) 进行 检测 。 利 用 Oligotex mRNA Kits Midi 试剂 
盒 说 明 书 ,纯化 各 样品 总 RNA 中 的 mRNA。 以 10 
hg mRNA 作为 模板 ,Gsul-oligo dT 作为 反 转录 引物 ， 
用 1 000 U Superscript 了 Reverse Transcriptase 在 


42C FRA 1 h 合成 cDNA 第 1 链 ;随后 利用 高 碘 酸 
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钠 氧化 mRNA 的 5' 端 帽子 结构 ,并 连接 生物 素 ; 通 
过 Dynal M280 磁 珠 躲 选 连接 了 生物 素 的 mRNA/ 
cDNA ,并 通过 碱 裂 解释 放 cDNA 第 1 链 ;然后 通过 
DNA Ligase 在 cDNA 第 1 链 的 5 末端 加 上 接头 , 利 























模式 。 
1.8 Gene Ontology ( GO) 5 KEGG 代谢 通路 富 
集 分 析 

利用 WEGO 软件 (Ye et al., 2006) 对 上 调 表达 








用 Ex Taq polymerase 合成 cDNA 第 2 链 。 最 后 , 通 
过 Gsu 工 酶 切除 polyA 和 5 ' 端 接头 。 利 用 Ampure 
beads Xj E XR cDNA 进行 纯化 , cDNA 文库 通过 
TruSeq"" DNA Sample Prep Kit-Set A 进行 构建 和 
TruSeq PE Cluster Kit 进行 扩 增 。 上 述 9 个 样品 的 次 
度 测 序 委托 广州 基 迪 奥 生 物 科 技 有 限 公 司 完 成 , 测 
序 平台 为 lumina Hiseq 2500 ,采用 双 端 (PE125 ) 测 
序 。 本 研究 测 得 的 转录 组 数据 已 上 传 美国 国家 生物 
技术 信息 中 心 (NCBI) SRA 数据 库 , SRA 号 
为 SRA454366。 
1.6 测序 数据 的 质 控 与 定位 

对 于 下 机 数据 ,利用 Pen 脚本 去 除 含有 
adaptor 未 知 核 苷 酸 比例 大 于 5$% 和 低 质量 reads , 获 
得 有 效 读 段 (clean reads) 。 使 用 短 reads 比 对 工具 
bowtie( Langmead et al., 2009) 分 别 将 各 样品 的 clean 
reads Bi] ( mapping) 到 核糖 体 数据 库 ( 最 多 允许 5 
个 错 配 ) (https://www. ncbi. nlm. nih. gov/genome/? 
term = apis% 20mellifera ) , 去 除 映射 上 核糖 体 的 
reads ,将 保留 下 来 的 数据 进一步 映射 意 蜂 的 参考 基 
因 组 ( Amel, 4. 5, http://www. ncbi. nlm. nih. gov/ 
genome/487? genome _ assembly id = 22683), 利 用 
SOAP aligner/soap2 软件 ( Hurgobin, 2016) 将 未 映射 
上 意 蜂 基因 组 的 reads 进而 映射 到 前 期 已 组 装 的 球 
SETS He RAOK, 2017 ) 。 
1.7 胁迫 意 蜂 幼虫 肠 道 的 球 圳 菌 DEGs 的 趋势 
分 析 

利用 FPKM (fragments per kilobase of transcript 
per million mapped reads) 法 计算 基因 表达 量 。 利 用 
R 软件 (version 2. 16. 2) 计算 各 样品 之 间 的 相关 性 
系数 。 利 用 edgeR 软件 (Robinson et al., 2010) 进行 
DEG 分 析 ( AamT2 vs AamT1 和 AamT3 vs AamT1) 。 
DEG HILERRIEN FDRs0.05 H.llog; Fold changel 
zl, 

为 了 对 各 样品 4, 5 和 6 日 龄 3 个 时 间 点 的 基 
因 表 达 模 式 进行 聚 类 ,首先 将 基因 表达 量 数据 归 一 
化 为 log,( 4AamT1/AamT1) , log, ( AamT2/ AamT1 ) 和 
log, ( AamT3/ AamTI1 ) ,利用 STEM (Short Time-series 
Expression Miner, vl. 3. 8) 软件 ( Ernst and Bar- 
Joseph, 2006) 进一步 分 析 基 因 表 达 模 式 。 在 所 有 的 
基因 表达 模式 中 ,P-value <0. 05 的 定义 为 显著 表达 















































模式 和 下 调 表 达 模 式 中 的 差异 表达 基因 ( DEGs) 进 
行 GO 富 集 分 析 。 利 用 Blastall 将 上 调 表达 模式 和 
下 调 表 达 模 式 中 的 DEGs 比 对 KEGG ( Kyoto 
Encyclopaedia of Genes and Genomes ) 数据 库 , 进行 
KEGG 代谢 通路 富 集 分 析 。 
1.9  RNA-seq 数据 的 实时 菊 光 定量 PCR ( RT- 
qPCR ) 验证 

为 了 验证 RNA-seq 数据 , 随机 选取 6 个 DEGs 
进行 RT-qPCR。 反 应 按照 SYBR Green Dye 试剂 盒 
(Vazyme xu], P Es) 操作 说 明 书 进行 。 本 实验 进 
行 3 次 生物 学 重复 和 3 次 技术 重复 。20 uL 的 反应 
体系 : SYBR Green Dye 10 uL,10. 0 umol /L E, 
向 引物 各 1 uL, cDNA 模板 DNA 1 pL, DEPC 水 7 
pL, RT-qPCR 反应 在 ABI 7500 Xt Œ PCR 4X 
(Life 公司 , 美国) 上 进行 ,反应 条 件 :95Y 预 变 性 1 
min;95Y 变性 15 s,60Y 延伸 30 s, 共 40 个 循环 ;最 
后 72% 延伸 45 s。 所 选 基因 的 相对 表达 量 采 用 
2 ( Livak and Schmittgen, 2001) 计 算 。 











2 结果 


2.1  RNA-seq 数据 质 控 与 评估 

球 昨 菌 胁迫 意 蜂 幼虫 肠 道 样品 的 Humina 测序 
共 得 到 球赛 菌 2$ 454 076 条 raw reads, 进一步 去 除 
低 质 数据 得 到 24 909 820 条 clean reads, 各 样品 
clean reads 数 均 在 2 253 412 (96.0896) 以 上 ,两 端 
Q20 均 在 97.08% 以 上 ,两 端 Q30 均 在 93. 46% 以 上 
( 表 1) ,上 述 结果 说 明 本 研究 测 得 的 转录 组 数据 质 
量 良 好 ,可 用 于 进一步 分 析 。 

进一步 对 AamT1，AamT2 和 AamT3 进行 主 成 
分 分 析 (PCA) ,结果 显示 第 一 主 成 分 (PC1 ) 与 第 二 
主 成 分 (PC2 ) 共 能 分 别 解 释 样品 基因 表达 总 体 方差 
的 88.2% 和 5.8% (图 1), 各 样品 的 组 内 重复 性 较 
好 ,而 不 同样 品 间 的 基因 表达 模式 总 体 差 异 明显 。 

RNA-seq 数据 比 对 球 喜 菌 参 考 转录 组 情况 统计 
显示 ,各 样品 clean reads 比 对 上 参考 转录 组 的 比例 
均 在 88.5296 以 上 ,AamT1 ，AamT2 与 AamT3 的 表 
达 基 因数 目 分 别 为 13 881, 17 682 和 20 026 个 ,分 
别 占 球 喜 菌 参考 转录 组 的 37.28% ，41.50% 和 
47.0196 (3:2), 

























































































apis-stressed 4- , 5- and 6-day-old larval gut samples of A. m. ligustica , respectively ( each with three replications, only 3-day-old larvae were fed with the 





artificial diet containing A. apis). 下 同 The same below. 
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E 肠 道 各 样品 的 球 溃 菌 的 PCA 分 析 
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Fig. 1 PCA analysis of Ascosphaera apis stressing larval gut samples of Apis mellifera ligustica 


4$ 陈 大 福 等 : 胁迫 意大利 蜜蜂 幼虫 肠 道 的 球赛 菌 的 转录 组 分 析 
表 1 胁迫 意 蜂 不 同日 龄 幼虫 肠 道 的 球 圳 菌 的 RNA-seq 数据 统计 
Table 1 RNA-seq datasets of Ascosphaera apis stressing different instar larval gut samples of Apis mellifera ligustica 
样品 。 原始 读 段 。 有效 读 段 ( 占 原始 读 段 的 百分比 ) QUU AEEA] UMORE 
Samples se Clear reads (percent uk rir ead umber and proportion of the Number and proportion of the 
99% correct bases ( Q20) 99.9% correct bases ( Q30) 
AamTl-1 2 528 840 2 469 622 (97.6696) 300 833 677 (97.4596) 291 359 297 (94.3896) 
AamT1-2 2 737 168 2 693 106 (98.3996) 329 389 184 (97.8596) 319 906 239 (95.0396) 
AamT1-3 3 023 132 2 980 240 (98.5896) 365 758 990 (98.1896) 356 656 737 (95.7496) 
AamT2-1 3 178 228 3 122 140 (98.2496 ) 382 437 729 (97.9996) 372 349 837 (95.4196) 
AamT2-2 2 547 204 2 501 890 (98.2296 ) 305 698 171 (97.7596) 296 488 524 (94.8096 ) 
AamT2-3 2 814 364 2 768 034 (98.3596) 339 309 805 (98.0796 ) 330 557 642 (95.5496) 
AamT3-1 3 229 530 3 178 978 (98.4396) 388 999 973 (97.8996) 378 044 002 (95.1496) 
AamT3-2 2 345 336 2 253 412 (96.0896 ) 273 461 324 (97.0896) 263 268 347 (93.4696) 
AamT3-3 3 050 274 2 942 398 (96.4696) 357 479 492 (97.1996) 344 314 704 (93.6196) 
AamTl, AamT2, AamT3; 分 别 是 球赛 菌 胁迫 的 意 蜂 4, 5 和 6 日 龄 幼虫 肠 道 样品 (各 3 个 重复 , 仅 3 日 龄 幼虫 饲 喂 仿 球赛 菌 的 人 工 饲料 ) A. 
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表 2 胁迫 意 蜂 幼虫 肠 道 各 样品 的 球 圳 菌 的 
表达 基因 数量 统计 
Table 2 Summary of the expressed genes in 
Ascosphaera apis stressing larval 


gut samples of Apis mellifera ligustica 





样品 表达 基因 数 表达 基因 比例 (% ) 
Samples ^ Number of expressed genes Proportion of expressed genes 
AamTI 15 881 37.28 

AamT2 17 682 41.50 

AamT3 20 026 47.01 


A Profile4: 0 P-value Profile7: 0 P-value 








Profile2: 1 P-value Profile3: 1 P-value 


Profile6: 0.0038 P-value 





2.2 ”胁迫 意 蜂 幼虫 肠 道 的 球 圳 菌 DEGs 的 趋势 分 析 

利用 STEM 软件 对 AamTl1，AamT2 和 AamT3 
的 表达 基因 进行 趋势 分 析 , 结果 显示 19 893 个 
DEGs 共聚 类 为 8 个 表达 模式 (profile0 ，profilel ， 
profile2, ， profile3 ， profile4 ， profileó 和 
profile7 ) (图 2) ,其 中 有 12 151 个 DEGs 聚 类 为 3 个 
表现 显著 上 升 趋势 的 表达 模式 ( profile4 ，profile6 和 
profile7) ( P-value <0.05)( 图 2)。profile4 包含 5 
953 个 DEGs, profile6 包含 2557 个 DEGs ,profile7 包 
含 3 641 个 DEGs。 





profile5 , 





Profile0: 0.07 P-value Profilel: 1 P-value 
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图 2 35 Eder dug e ESTERI. DEG 的 趋势 分 析 
Fig. 2 Trend analysis of DEGs in Ascosphaera apis stressing larval gut samples of Apis mellifera ligustica 
A: 趋势 分 析 结 果 Trend analysis result; B; profile 包含 DEGs 的 表达 模式 Expression patterns of profile4; C; profile6 包含 DEGs 的 表达 模式 
Expression patterns of profile6 ; D: profile7 包含 DEGs 的 表达 模式 Expression patterns of profile7. 


2.3 GO BENIT 

上 调 趋 势 中 所 有 DEGs 的 GO 分 类 结果 显示 ,这 
些 DEGs 分 为 3 类 :生物 学 进程 (biological process ) , 
细胞 组 分 (cellular component ) 和 分 子 功 能 (molecular 
function) ,它们 富 集 于 40 个 GO term ,分 别 为 细胞 进 
程 (cellular process) 代谢 进程 (metabolic process) 、 细 
胞 (cell) 、 细 胞 组 件 (cell part ) , ff& f E ( catalytic 
activity) 单 组织 进 程 (single-organism process) 结合 
(binding )、 细 胞 器 (organelle )、 大 分 子 复合 物 
( macromolecular complex )、 细 胞 器 组 件 ( organelle 
part), Æ # J 控 ( biological regulation )、 定 位 


(localization ) 、 细 胞 膜 (membrane ) 、 细 胞 成 分 组 织 
(cellular component organization or biogenesis) 、 细 胞 膜 
组 件 (membrane part) 、 应 激 (response to stimulus) 、 转 
运 器 活性 (transporter activity)、 结 构 分 子 活 性 
(structural molecule activity ) , Jl 内 f$ ( membrane- 
enclosed lumen )、 信 * (signaling), / 8 ut fi 
( developmental process ) , ^E 7É ( reproduction ) , 生殖 进 
fi ( reproductive process) | Tz Tr 2 28 cr SK PS T i TE 
(nucleic acid binding transcription factor activity ) , 4) 
子 功 能 调节 器 (molecular function regulator ) , ^: K 
(growth) 多 组 织 进程 (multi-organism process ) | Jj BE 
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粒子 (virion) ,病毒 粒子 组 件 (virion part) 、 多 细胞 生 
物 进程 (multicellular organismal process) 分 子 传感器 
活性 (molecular transducer. activity )、 抗 氧化 活性 


号 传感器 活性 (signal transducer activity ) .解毒 作用 
( detoxification ) 、 胞 外 区 (extracellular region ) 细胞 连 
接 (cell junction )、 生物 附 着 (biological adhesion ) 











(antioxidant activity )、 转 录 因 子 结合 蛋白 活性 
(transcription factor activity, protein binding) .电子 载 
体 活 性 (electron carrier activity)、 拟 核 (nucleoid) 、 信 


(图 3)。 富 集 基 因数 最 多 的 CO term 为 细胞 进程 
(2 601 unigenes) ,其 次 为 代谢 进程 (2 553 unigenes ) 
和 细胞 (2 522 unigenes) (图 3)。 
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图 3 胁迫 意 蜂 幼虫 肠 

















道 各 样品 的 球 宫 菌 DEGs 中 表达 上 调 基因 的 GO 富 集 分 析 


Fig. 3 GO enrichment analysis for the up-regulated expressed DEGs in Ascosphaera apis stressing 

larval gut samples of Apis mellifera ligustica 
1: 生物 附着 Biological adhesion; 2: 生物 调控 Biological regulation; 3: 细胞 组 分 组 织 或 
细胞 进程 Cellular process; 5: 解毒 作用 Detoxification; 6: 发 育 进 程 Developmental process; 7: 生长 Growth; 8: 定位 Localization; 9: 代谢 进程 
Metabolic process; 10; 多 组 织 进程 Multi-organism process; 11: 多 细胞 生物 进程 Multicellular organismal process; 12; 生殖 Reproduction; 13: 生殖 
进程 Reproductive process; 14; 应 激 Response to stimulus; 15; 信号 Signaling; 16; 单 组 织 进程 Single-organism process; 17: 细胞 Cell; 18 : 细胞 连 
接 Cell junction; 19; 细胞 组 件 Cell part; 20; 胞 外 区 Extracellular region; 21: 大 分 子 复 合 物 Macromolecular complex; 22: 细胞 膜 Membrane; 23 : 





生物 合成 Cellular component organization or biogenesis; 4; 
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8: 病毒 粒子 Virion; 29; 





























Structural molecule activity; 39; 





2.4 KEGG 代谢 通路 富 集 分 析 

对 上 调 趋势 中 所 有 DEGs 进行 KEGG 代谢 通路 
富 集 分 析 ,结果 显示 这 些 DEGs 富 集 在 119 个 代谢 
通路 上 , 其 中 富 集 基因 数量 最 多 的 是 核糖 体 
(ribosome) (213 unigenes) ,其 次 为 氨基 酸 生 物 合成 
( biosynthesis of amino acids) (154 unigenes ) 和 内 质 
网 蛋白 加 工 (protein processing in endoplasmic 
reticulum) (130 unigenes) (图 4, 表 3)。 富 集 基 因 
数 前 20 位 的 代谢 通路 如 表 3 所 示 , 有 大 量 DEGs 富 
集 在 各 类 物质 代谢 , 如 碳 代谢 (carbon metabolism ) 
(129 unigenes)、 味 叭 代谢 (purine metabolism ) ( 119 
unigenes)、 喀 啶 代谢 (pyrimidine metabolism ) ( 86 
unigenes) 以 及 半 胱 氨 酸 与 蛋氨酸 代谢 (cysteine and 

















胞 膜 组 件 Membrane part; 24: 细胞 膜 内 腔 Membrane enclosed lumen; 25; 拟 核 Nucleoid; 26 : 
病毒 粒子 组 件 Virion part; 30: 抗 氧 化 活性 Antioxidant activity; 31 : 

电子 载体 活性 Electron carrier activity; 34: 分 子 功 能 调节 因子 Molecular function regulator; 35 : 
核酸 结合 转录 因子 活性 Nucleic acid binding transcription factor activity; 37: 信和 
转录 因子 活性 , 蛋白 结合 Transcription factor activity, protein binding; 40; 转运 活性 Transporter activity. 


细胞 器 Organelle; 27: 细胞 器 组 件 Organelle part; 
结合 Binding; 32: 催化 活性 Catalytic activity; 33; 
分 子 传感器 活性 Molecular transducer activity ; 36 : 


号 传感器 活性 Signal transducer activity; 38: 结构 分 子 活 性 


methionine metabolism) (57 unigenes) 。 共 有 48 个 
DEGs 富 集 在 MAPK 信和 号 通路 上 ,对 富 集 在 MAPK 
言 号 通路 上 的 基因 进行 表达 量 聚 类 ,结果 显示 随 着 
胁迫 时 间 的 延长 , 富 集 在 该 通路 上 的 基因 表达 水 平 








显著 增强 (图 5) 。 
2.5  RNA-seq Zi ja rh Ek X€ f] DEGs 表达 水 平 的 
RT-qPCR 验证 


利用 RT-qPCR 检测 随机 选取 的 6 个 DEGs, £1 
括 1 个 组 蛋白 HA. 1 编码 基因 、1 个 Pe20g13990 4£& 
白 编码 基因 、3 个 未 知 基因 和 1 个 醛 酮 还 原 酶 编码 
基因 。 结 果 显 示 这 些 DEGs 的 表达 水 平 变化 趋势 与 
RNA-seq 数据 中 的 基因 表达 水 平 变化 趋势 一 致 (图 
6) ,证 明 本 研究 获得 的 转录 组 数据 真实 可 乱 。 
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图 4 ”胁迫 意 蜂 幼虫 肠 道 各 样品 的 球 填 菌 DEGs 表达 上 调 基 因 的 KEGG 代谢 通路 富 集 分 析 
Fig. 4 KEGG pathway enrichment analysis for the up-regulated expressed DEGs in Ascosphaera apis 
stressing larval gut samples of Apis mellifera ligustica 
表 3 上 升 趋势 中 富 集 基因 数 前 20 位 的 代谢 通路 
Table3 Top 20 pathways enriched by unigenes within the up-regulated expression profiles 
代谢 通路 基因 数 代谢 通路 TD 
. P-value Q-value 
Metabolic pathway Number of genes Pathway ID 

核糖 体 Ribosome 213 5.09E-06 0.000303 ko03010 
氨基 酸 生 物 合成 Biosynthesis of amino acids 154 0.166726 0.901836 ko01230 
内 质 网 蛋白 加 工 Protein processing in endoplasmic reticulum 130 0.614214 0.997893 ko04141 
碳 代 谢 Carbon metabolism 129 0.921245 0.997893 ko01200 
味 叭 代谢 Purine metabolism 119 0. 522226 0.997893 ko00230 
剪接 体 Spliceosome 119 0.847405 0.997893 ko03040 
细胞 周期 -酵母 Cell cycle - yeast 114 0.000156 0.006176 ko04111 
减 数 分裂 - 酵 母 Meiosis - yeast 111 2.53E-06 0.000301 ko04113 
氧化 磷酸 化 Oxidative phosphorylation 111 0.369017 0.997893 ko00190 
RNA 转运 RNA transport 111 0.899856 0.997893 ko03013 
内 春 作 用 Endocytosis 103 0.818927 0.997893 ko04144 
泛 素 介 导 的 蛋白 水 解 Ubiquitin mediated proteolysis 91 0.414170 0.997893 ko04120 
喀 喧 代谢 Pyrimidine metabolism 86 0.559692 0.997893 ko00240 
真 核 生物 核糖 体 生 物 合成 Ribosome biogenesis in eukaryotes 85 0.626941 0.997893 ko03008 
RNA 降解 RNA degradation 72 0.539455 0.997893 ko03018 
糖 醇 解 / 糖 异 生 Glycolysis/Gluconeogenesis 69 0.263767 0.997893 ko00010 
mRNA 监视 代谢 通路 mRNA surveillance pathway 65 0.432491 0.997893 ko03015 
TW P Phagosome 63 0.772199 0.997893 ko04145 
核 苷 酸 切 除 修复 Nucleotide excision repair 57. 0. 129344 0.810104 ko03420 
半 胱 氨 酸 与 蛋氨酸 代谢 Cysteine and methionine metabolism 57 0. 19747 0.903804 ko00270 
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图 5 MAPK 信号 通路 富 集 基因 的 热 图 
Fig. 5 Heatmap of unigenes enriched in MAPK signaling pathway 
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6 胁迫 意 蜂 幼虫 肠 道 各 样品 的 球赛 菌 RNA-seq 数据 的 RT-qPCR 验证 
Fig. 6 RT-qPCR validation of RNA-seq data of Ascosphaera apis stressing larval gut samples of Apis mellifera ligustica 
A: 组 蛋白 H4. 1 编码 基因 Histone H4. 1 encoding gene; B; Pc20g13990 和 蛋白 编码 基因 Pc20813990 encoding gene; C - E: 未 知 基因 Unknown 
gene; F: je ddp] Aldo-keto reductase encoding gene. 

















参考 转录 组 ( 张 盟 楠 等 , 2017 ) ,为 球 喜 菌 的 转录 组 
3 Wit 学 研究 英 定 了 基础 。 在 Comman 45 (2012) 的 研究 
中 ,测序 对 象 是 来 自 培养 基 的 球 圳 菌 菌 丝 和 来 自 密 

此 前 ,本 课题 组 de novo 组 装 并 注释 了 球 圳 菌 的 。” 蜂 幼虫 感染 组 织 的 球 圳 菌 菌 丝 , 其 中 后 者 是 感染 幼 
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虫 肠 道 内 萌发 的 菌 丝 穿 透 肠 壁 后 继续 生长 最 后 突破 
体 表 的 苗 丝 ,其 转录 组 变化 并 不 能 精确 反映 侵 染 过 
程 中 球 宫 菌 的 基因 表达 情况 ,因为 此 时 的 菌 丝 已 经 
不 与 宿主 互 作 。 本 研究 的 测序 对 象 胁迫 意 蜂 幼虫 肠 
道 的 球 澡 菌 , 肠 道内 的 球 吉 菌 通过 与 宿主 的 互 作 促 
进 侵 染 过 程 , 故 通过 映射 意 蜂 参考 基因 组 过 滤 得 到 
的 球 吉 菌 转录 组 数据 ,能 更 精确 地 反映 病原 在 侵 染 
过 程 中 的 基因 表达 谱 。 

研究 中 ,趋势 分 析 结 果 显示 ,胁迫 的 意 蜂 幼虫 肠 道 的 
球 宫 菌 DEGs 的 显著 表达 模式 丝 为 上 调 趋 势 ,说 明 
在 胁迫 意 蜂 幼虫 肠 道 的 过 程 中 ,球赛 菌 的 基因 在 很 
大 程度 上 被 激活 。 本 研究 发 现在 胁迫 意 蜂 幼虫 肠 道 
过 程 中 , 球 吉 苗 参与 物质 代谢 的 多 数 基 因 均 表现 出 
上 调 趋势 ,大 量 DEGs 富 集 在 多 达 80 个 物质 代谢 通 






















































































内 质 网 蛋白 质 加 工 代 谢 通路 中 的 泛 素 家 族 重 白 编码 
基因 (unigene0001921 ) ATP 4K ift KI ^ F ££ 4H 
HSC82 ( unigene0002086 ) 及 热 激 和 蛋白 90-2 编码 基 
(unigene0036009 ) 等 。 

真菌 可 通过 精细 而 保守 的 信号 通路 ,如 调节 黑 
色素 形成 .配合 、 毒 力 以 及 形态 发 生 等 来 感知 外 界 环 
境 并 迅速 做 出 应 答 ( Cornman et al., 2012)。 在 真菌 
的 信号 通路 中 ,MAPK 信和 号 通路 最 为 关键 ,在 真菌 致 
病 性 及 响应 环境 胁迫 方面 发 挥 重要 作用 (Perez and 
Cansado, 2010)。 本 研究 发 现 48 个 DEGs 富 集 在 
MAPK 信号 通路 上 , 如 MAPK MKCI 编码 基 
(unigene0001632 ) , MAPK 激酶 MKK1/SSP32 编码 
基因 (unigene0003669 ) 以 及 核 孔 复合 物 蛋 白 An- 
Mlpl 编码 基因 (unigene0039803 ) 等 。 本 研究 还 发 现 
EX DEGs 的 表达 水 平 随 着 球 宫 菌 胁迫 时 间 的 延长 










































































路 ,包括 碳水 化 合 物 代 谢 、 脂 代谢 和 氨基 酸 代谢 等 。 
伴随 物质 代谢 的 增强 , 球 吉 菌 的 能 量 代 谢 在 侵 染 过 











而 显著 增强 (图 5 ) ,推测 该 通路 在 球赛 菌 对 意 蜂 幼 
虫 的 致 病 性 方面 发 挥 关 键 作 用 ,通过 RNAi 技术 对 

















程 中 显著 增强 ,表现 为 参与 能 量 代 谢 的 多 数 基因 旦 
现 上 调 趋势 ,这 些 基因 富 集 在 4 个 能 量 代谢 通路 上 ， 
如 氧化 磷酸 化 代谢 通路 中 的 V 型 ATP 酶 了 上 亚 基 编 
码 基 因 (unigene0000213 ) , NADH 脱氧 酶 亚 基 5 编 
码 基因 (unigene0002084) 和 NADH- 辅 酶 Q 氧化 还 
原 酶 编码 基因 (unigene0003823 ) SE, PEW ER ZE P AN 
意 蜂 幼虫 通过 病原 -宿主 互 作 使 二 者 相互 适应 , 球 吉 
菌 在 宿主 肠 道内 通过 增强 物质 和 能 量 代谢 来 促进 侵 
染 过 程 。 

在 胁迫 意 蜂 幼虫 肠 道 的 过 程 中 , 球 法 菌 需 不 断 
合成 遗传 物质 和 和 蛋白质, 从 而 满足 增殖 需要 ,本 研究 
发 现 大 量 参与 遗传 物质 与 蛋白 质 合 成 的 基因 显 若 上 
调 , 如 DNA 复制 (DNA replication ) 中 的 复制 因子 C 
亚 基 1 编码 基因 (unigene0002253) .DNA 复制 许可 
因子 MCM4 编码 基因 ( unigene0004671 ) 及 核糖 核酸 
Hg H2 亚 基 C 编码 基因 (unigene0006412) ; RNA 转 
i5 ( RNA transport) 中 的 翻译 起 始 因子 3 亚 基 丁 编码 
基因 (unigene0002179 ) 核糖 核酸 酶 复合 物 亚 基 
Pop4 编码 基因 ( unigene0004409 ) 及 精 氨 酸 N- 甲 基 
转移 酶 Skbl 编码 基因 (unigene0006216 ) 等 ;核糖 体 
中 的 eos 核糖 体 蛋 白 LU 编码 基因 
( unigene0000035 ) .60S 核糖 体 蛋白 L21-A 编码 基因 
(unigene0000859 ) 及 60S 核糖 体 蛋 白 LIS 编码 基因 
( unigene0041982 ) 等 ;氨基 酸 生 物 合成 代谢 通路 中 
的 丙酮 酸 羧 化 酶 编码 基因 (unigene0024192 ) ,3-875 
酸 甘油 醛 脱 氧 酶 编码 基因 (unigene0008207 ) 及 果 
糖 -二 磷酸 醛 缩 酶 A 编码 基因 (unigene0041009 ) 等 ; 





















































































































































该 通路 上 的 基因 进行 沉默 ,或 许 是 治疗 白垩 病 的 一 
种 可 行 途径 。 

本 人 研究 利用 RNA-seq 技术 对 球 圳 菌 胁迫 的 意 
蜂 幼 虫 肠 道 进 行 深度 测序 ,基于 趋势 分 析 进 一 步 对 
胁迫 过 程 中 的 球 喜 菌 进 行 转录 组 学 分 析 ,研究 结果 
为 在 分 子 水 平 揭示 球赛 苦 的 致 病 机 理 提 供 了 重要 信 
息 ,也 为 阐释 白垩 病 过 程 中 的 病原 -宿主 互 作 机 制 黄 
定 了 基础 。 
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